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ABSTRACT 

Context. The gas-driven dust activity of comets is still an unsolved question in cometary physics. Homogeneous dust layers com¬ 
posed of micrometer-sized grains possess tensile strengths of ~ 1 kPa, which is far higher than typical gas pressures caused by the 
sublimation of the ices beneath the covering dust layer. This implies that the dust grains cannot be detached from the surface by 
the gas pressure of the sublimating ices. One possibility to avoid this problem is that the nucleus formed through the gravitational 
collapse of an ensemble of millimeter- to centimeter-sized aggregates. In this case, an aggregate layer with a tensile strength on the 
order of ~ 1 Pa is formed on the surface of the nucleus, which allows for the release of the aggregates from the surface by the gas 
pressure build up at the ice-dust interface. 

Aims. We use the gravitational instability formation scenario of cometesimals to derive the aggregate size that can be released by the 
gas pressure from the nucleus of comet 67P/Churyumov-Gerasimenko for different heliocentric distances and different volatile ices. 
Methods. To derive the ejected aggregate sizes, we developed a model based on the assumption that the entire heat absorbed by the 
surface is consumed by the sublimation process of one volatile species. The calculations were performed for the three most prominent 
volatile materials in comets, namely, H 2 O ice, CO 2 ice, and CO ice. 

Results. We find that the size range of the dust aggregates able to escape from the nucleus into space widens when the comet 
approaches the Sun and narrows with increasing heliocentric distance, because the tensile strength of the aggregates decreases with 
increasing aggregate size. The activity of CO ice in comparison to H 2 O ice is capable to detach aggregates smaller by approximately 
one order of magnitude from the surface. As a result of the higher sublimation rate of CO ice, larger aggregates are additionally able 
to escape from the gravity field of the nucleus. 

Conclusions. Our model can explain the large grains (ranging from 2 cm to 1 m in radius) in the inner coma of comet 67P/Churyumov- 
Gerasimenko that have been observed by the OSIRIS camera at heliocentric distances between 3.4 AU and 3.7 AU. Furthermore, the 
model predicts the release of decimeter-sized aggregates (trail particles) close to the heliocentric distance at which the gas-driven dust 
activity vanishes. However, the gas-driven dust activity cannot explain the presence of particles smaller than ~ 1mm in the coma 
because the high tensile strength required to detach these particles from the surface cannot be provided by evaporation of volatile ices. 
These smaller particles can he produced for instance by spin-up and centrifugal mass loss of ejected larger aggregates. 

Key words. Comets: general. Methods: analytical. Solid state: volatile 


1. Introduction 

Comets and their activity are most often associated with the 
emission of dust that forms a visible coma and a tail, the shape 
of the latter determined by the radiation pressure of the sun. The 
role of volatiles for cometary activity, first observed by the emis¬ 
sion of radicals but not of stable molecules, was initially strongly 
underestimated, which led to models such as the sandbank model 
(Lyttleton 1953). Better spectroscopic observations and the de¬ 
termination of nongravitational forces led to the dirty snowball 
model (Whipple 1950, 1951), where ices such as solid H 2 O, CO 2 
, and CO, polluted by dust, form a solid nucleus. The dust is then 
liberated by the sublimation of the abundant volatile ices. The 
first space missions to comet IP/Halley (Sagdeev et al. 1986; 
Keller et al. 1986) revealed a pitch-black nucleus with a density 
of only 0.6 )^q 4 g cm“^. It was determined that dust, that is, a re¬ 
fractory material, is more abundant than volatiles. In situ mea¬ 
surements and estimates revealed a dust-to-gas mass ratio some¬ 


where between 1 and 10. The observations of comet IP/Halley 
also showed that only a small part of its surface was active. Most 
of the surface seemed to be covered by an inert crust for which 
cohesive forces between the dust grains are much stronger than 
gravitational forces (Heim et al. 1999). The question how activ¬ 
ity works and how it is maintained over many cometary orbits 
remained unanswered. 

Modeling efforts concentrated on the heat transport into the 
porous nucleus, which is required to sublime the ices. The ma¬ 
terial structure (porosity and packing structure) and the physi¬ 
cal properties (e.g., heat conductivity) strongly influence the gas 
production rates (Gundlach et al. 2011). Modeling parameters 
were adjusted to fit the observations (see, e.g., Prialnik et al. 
2004; Enzian 1999). In the models, the dust was somehow lib¬ 
erated by the gas, leaving the surface, and was accelerated in 
the inner coma against the minute gravity of the nucleus (Ein- 
son & Probstein 1968). Kuehrt & Keller (1994) showed that the 
sublimation vapor pressure is much too weak to overcome the 
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van-der-Waals forces between small dust grains. Following this 
paper, Mohlmann (1995) presented typical values of cohesion 
of porous media and concluded that a cohesively bound matrix 
of relatively refractory and porous matter cannot be destroyed by 
the gas drag. To work around this insuperable obstacle, the cohe¬ 
sion force between microscopic mineral grains was just ignored 
in most similar publications (see, e.g., Rickman et al. 1990). 
Once this dust crust has formed, activity is quenched by the re¬ 
duction of both the heat flow into the interior and the gas flux 
from the volatile subsurface region of the cometary nucleus. 

Recently, however, it was shown that comets can form by 
gravitational collapse of centimeter-sized dust agglomerates in 
the early solar system (Johansen et al. 2007; Skorov & Blum 
2012; Blum et al. 2014, 2015). There are strong indications that 
the dust and water-ice agglomerates were mixed by radial trans¬ 
port of the particles before the cometesimals formed in the outer 
solar system beyond the CO sublimation boundary (see results 
of the Stardust mission; Brownlee et al. 2006). The gravitational 
collapse then forms a cometesimal composed of porous dust and 
ice aggregates, themselves consisting of micrometer-sized parti¬ 
cles (Sunshine et al. 2007, 2011; Protopapa et al. 2014). The ex¬ 
istence of separate ice and dust aggregates composed of micron¬ 
sized particles is also suggested by the observations of grains 
containing water ice in the inner coma of comet 103P/Hartley 2 
(A’Hearn et al. 2011) and on the surface of comet 9P/Tempel 1 
(Sunshine et al. 2006). 

When the comet approaches the Sun for the first time, the 
ice on the surface sublimes and the remaining dust agglomerates 
form a thin, non-volatile dust layer (Skorov & Blum 2012). The 
ice-free surface material possesses an ultra-low tensile strength 
(i.e., the cohesive force per unit area) of only a few Pascal, or 
less in the case of millimeter-sized aggregates (Blum et al. 2014). 
This is similar to the highest gas pressures that can be reached 
by the outgassing of typical volatile constituents such as H 2 O 
ice, CO 2 ice, or CO ice, at typical subsurface temperatures of 
comet nuclei (Blum et al. 2014, 2015). The thin layer of dust 
aggregates can then be ejected if the gas pressure beneath it ex¬ 
ceeds the tensile strength of the material. The vapor pressure 
of H 2 O ice suffices to support the dust activity of such a gen¬ 
tly formed nucleus at heliocentric distances smaller than 2.5 AU 
(Blum et al. 2015). The aggregates are released intact because of 
their relatively high internal tensile strength (Blum et al. 2006; 
Kothe et al. 2010). Icy aggregates will have to be lifted with 
the help of super-volatiles (e.g., CO 2 ice, or CO ice) in agree¬ 
ment with the above-mentioned observations in the inner coma 
of comet Hartley 2. Blum et al. (2014) have already shown that 
the ejection of H 2 O ice aggregates is formally possible by the 
outgassing of H 2 O ice at 1 AU (see their Fig. 10). However, 
as this is physically not possible, materials possessing higher 
volatility are required to release H 2 O ice aggregates. 

At larger heliocentric distances, the sublimation pressure of 
H 2 O ice is too low to explain the dust emission of comets, how¬ 
ever. Beyond 2.5 AU, more volatile icy compounds are needed 
to explain the observed activity of comets. CO 2 and in particular 
CO ice possesses a lower latent heat (i.e., a higher sublimation 
rate if the same amount of energy is available for sublimation) 
than H 2 O. Indeed, CO 2 and CO are frequently observed in rela¬ 
tively high abundances of ~ 10 % to ~ 20 % with respect to H 2 O 
(see, e.g., Bockelee-Morvan et al. 2004; A’Heam et al. 2012). 
At these heliocentric distances, H 2 O ice is inactive, so that the 
outgassing of CO 2 ice, or CO ice causes the ejection of H 2 O ice 
aggregates. The ejection of H 2 O ice aggregates by the activity of 
CO 2 , or CO ice is also possible at shorter heliocentric distances. 


however, as can be seen in the case of comet 103P/Hartley 2 
(AHearnet al. 2011). 

The observed CO production rate of comet Hale-Bopp 
(C/1995 Ol) directly followed the variation of the solar irradi¬ 
ation from 5 AU down to perihelion and out again (Biver et al. 
1999; Gortsas et al. 2011). This suggests that CO ice must be 
present relatively close to the surface of the nucleus even at per¬ 
ihelion. In the case of comet 103P/Hartley 2, the outgassing of 
CO 2 ice was observed to drive the ejection of centimeter-sized 
icy particles (A’Hearn et al. 2011). These observations strongly 
support the hypothesis that super-volatiles can survive close to 
the eroding surface and might be real drivers of observable dust 
activity. 

To fit numerous photometric and colorimetric observations 
of comet dust (Kolokolova et al. 2007), cometary grains have 
to be porous aggregates. Comet polarimetry in particular has 
shown that these aggregates in turn consist of micron and sub¬ 
micron porous pieces. This structure of the cometary particles 
leads to a dramatic reduction of the van-der-Waals forces be¬ 
tween the contacting aggregates (Skorov & Blum 2012). As 
noted previously (Skorov & Blum 2012; Blum et al. 2014), it 
is easier to lift dust layers consisting of large agglomerates than 
those fromed by small ones. Thus the evaluation of cohesion 
between the aggregates provides the lower limit of the size of 
the ejected aggregates. On the other hand, the gravity of the 
comet restricts the maximum size of the aggregates that can be 
accelerated away by the out-flowing gas (Whipple 1950). The 
early observations in the inner coma of comet 67P/Churyumov- 
Gerasimenko by the Rosetta spacecraft revealed a considerable 
amount of large particles. Images taken by the OSIRIS cam¬ 
era (Keller et al. 2007) even at heliocentric distances larger 
than 3.6 AU are full of tracks of particles of up to 2 cm in di¬ 
ameter (Rotundi et al. 2015). Furthermore, particles with di¬ 
ameters between 4 cm and 2 m were observed in bound orbits 
close to the nucleus of comet 67P/Churyumov-Gerasimenko be¬ 
tween 3.4 AU and 3.7 AU. These bound particles were possi¬ 
bly emitted during or after the last perihelion passage of comet 
67P/Churyumov-Gerasimenko. 

In this work, we discuss how the gas-driven dust activity of 
comets can be modeled based on laboratory results and on the 
assumption that comets have formed by gravitational instability 
(see Sect. 2). In Sect. 3, this model is then applied to comet 
67P/Churyumov-Gerasimenko to show that the outgassing of 
CO 2 ice and CO ice preferentially ejects large aggregates, which 
has been observed by Rosetta, and that evaporating H 2 O ice can¬ 
not eject any particles at large heliocentric distances (> 2.5 AU). 
Section 4 summarizes the main results of this work. 

2. Modeling cometary dust activity 

To release dust particles from a cometary nucleus, the force ex¬ 
erted on the dust by the outflowing gas molecules must exceed 
the cohesive force of the dust aggregates. Here, we recall that 
the assumption of an ice-free surface is an idealization of the 
problem, but observations have shown that most of the cometary 
surfaces are covered by a nonvolatile dust layer (see, e.g., Keller 
et al. 1988; Buratti et al. 2004). The tensile strength of the dust- 
aggregate layer reads 

y^Yq<ps-^i\ ( 1 ) 

where To = 1.6 Pa is a constant, (p is the volume filling factor of 
the aggregate layer, and s is the dust-aggregate radius measured 
in millimeters (Skorov & Blum 2012). Brisset et al. (2015) ex¬ 
perimentally showed that Eq. 1 is quantitatively correct (see their 
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Figure 1. Gas pressure beneath the covering dust layer as a function of the thickness of the dust layer (in units of the dust-aggregate diameter) 
at different heliocentric distances (between 1AU and 4 AU) and for different volatile materials (H 2 O ice, CO 2 ice and CO ice). The solid blue 
curves and the black dashed lines show the gas pressure derived by the thermophysical model and the approximation by Eq. 3, respectively. The 
temperatures of the ice surface at the peak of the pressure curves are 208.4 K (H 2 O ice at 1 AU), 196.3 K (H 2 O ice at 2 AU), 117.2 K (CO 2 ice at 
2AU), 106.6 K (CO 2 ice at 3 AU), 33.9 K (CO ice at 3 AU), and 32.8 K (CO ice at 4AU). The deviation of the approximation from the highest 
pressures derived by the thermophysical model are denoted by the arrows. The sublimation pressure is calculated by \n(ps„b) = a - bjT, with 
the coefficients a = 3.23 x lO'^ Pa and b = 6134.6 K taken from Gundlach et al. (2011) for H 2 O ice. For CO 2 ice and CO ice the coefficients are 
a = 1.23 X 10*^ Pa and b = 3167.8 K as well as a = 1.26 x 10® Pa and b = 764.2 K, respectively (Fanale & Salvail 1990). We used a Bond albedo 
of 0.04 for the calculations (the geometric albedo of comet 67P/Churyumov-Gerasimenko is 0.05; Sierks et al. 2015). 
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Fig. 17). The smallest dust-aggregate radius that can be released 
from the surface by the sublimation of the ices is given if the gas 
pressure p equals the tensile strength of the aggregate layer (Eq. 
1 ), 


^min 


p 

Y cpj 


( 2 ) 


For cometary surface layers formed by gravitational instabil¬ 
ity (i.e., for dust aggregate layers), the tensile strength is very 
low (~ 1 Pa) and somewhat counterintuitively decreases with in¬ 
creasing dust-aggregate size (Skorov & Blum 2012; Blum et al. 

2014) . Hence, it is easier to release larger dust aggregates than 
smaller ones. For example, a densely packed aggregate layer 
with a volume-filling factor of 0.6 consisting of lOO/rm-sized 
dust aggregates possesses a tensile strength of 4.5 Pa, whereas 
increasing the aggregate size to 1 cm yields a tensile strength of 
0.2 Pa. 

The gas pressure p of the sublimating ice, covered by a 
porous dust layer, depends on its transport characteristics (gas 
permeability) and on the energy supplied to the region. The lat¬ 
ter depends on the heliocentric distance, on the nucleus albedo, 
on the thickness, and on the heat conductivity of the covering 
dust layer. A detailed heat-conductivity model for granular ma¬ 
terials in vacuum was presented by Gundlach & Blum (2012) 
and Gundlach & Blum (2013). Together with the available en¬ 
ergy for the sublimation process, the permeability of the cover¬ 
ing dust layer determines the pressure build-up at the depth of 
the ice-dust interface (Gundlach et al. 2011; Blum et al. 2014, 

2015) . In the following, this detailed model is referred to as the 
thermophysical model. 

In general, the derivation of the gas pressure at the ice- 
dust interface should be determined from full consideration of 
the conservation of energy, momentum, and mass for solid and 
gaseous phases. However, such a detailed consideration is far 
beyond the scope of this paper. We wish to simplify the physical 
model as much as possible (keep it reasonable) to highlight the 
new qualitative results. For this reason, we simplified the deriva¬ 
tion of the gas pressure at the ice-dust interface by the approxi¬ 
mation that the incoming energy that is absorbed by the nucleus 
is entirely used for the sublimation of the volatiles (i.e., thermal 
reradiation of the comet nucleus and energy dissipation by its 
core are neglected). The reasonableness of this simplification is 
illustrated below. With this assumption, the gas pressure at the 
ice-dust interface is given by 


p -A)Ie 


j ink Tice 

fl/, / A V m 


(3) 


where A is the Bond albedo of the surface material, Ie is the 
solar constant, qe is the distance of the Earth to the Sun, a/, is 
the heliocentric distance of the comet, A is the latent heat of the 
sublimating ice, k is the Boltzmann constant, Tice is the tempera¬ 
ture of the evaporating ice, and m is the mass of the evaporating 
molecules. 

To illustrate the plausibility of this simplification, the results 
of Eq. 3 are compared with the corresponding pressure obtained 
from the sophisticated thermophysical model for H 2 O ice (latent 
heat of sublimation; A = 2.6 x 10® Jkg '; Orosei et al. 1995) 
at the ice-dust interface for different thicknesses of the covering 
dust-aggregate layer and at two different heliocentric distances 
(at 1 AU and 2 AU; see Figs, la and lb). 

For another test of the validity of the approximation, we cal¬ 
culated the resulting gas pressure for CO 2 ice (at 2 AU and 3 AU; 


latent heat of sublimation; A = 5.7 x 10® Jkg"^; Mavko et al. 
2009, see Figs. Ic and Id) and for CO ice (at 3 AU and 4 AU; la¬ 
tent heat of sublimation; A = 3.0 x 10® Jkg '; Brown & Ziegler 
1979, see Figs, le and If). These two volatile materials were 
chosen because of their relatively high abundance in comets and 
because model calculations suggest that volatile species incor¬ 
porated in planets after their formation are mainly composed of 
H 2 O, CO 2 , CO, CH 3 OH, and NH 3 (Marboeuf et al. 2014). 

A comparison of the derived gas pressure at the ice-dust in¬ 
terface using the thermophysical model (solid blue curves in Fig. 
1) with the approximation of Eq. 3 (dashed black line; derived 
for the temperature of the ice-dust interface at which the high¬ 
est pressure is reached in the thermophysical model) yields the 
deviations indicated by the arrows in Fig. 1. The pressure ra¬ 
tio between the approximation and the thermophysical model 
indicated in Fig. 1 was calculated for the maxima of the pres¬ 
sure curves, because the maximum yields the smallest aggregate 
size that can be released by the gas pressure. For the considered 
cases, the maximum gas pressure is typically reached for a dust 
layer thickness between two and five aggregate layers. 

Figure shows 1 that the deviation from the thermophysical 
model to the approximation of Eq. 3 increases with increasing 
heliocentric distance. This can be explained by the following 
considerations. For reasonable dust-aggregate sizes and surface 
temperatures, the heat conductivity is dominated by radiative 
transport and, thus, increases strongly with increasing tempera¬ 
ture, that is, with decreasing heliocentric distance (see Gundlach 
& Blum 2012; Gundlach & Blum 2013, for details). Therefore, 
the fraction of absorbed energy that is transported into the inte¬ 
rior of the comet nucleus and, thus, used to evaporate volatile 
ices, also increases with decreasing heliocentric distance. On 
the other hand, for increasing heliocentric distances, the relative 
importance of thermal reradiation of the surface into deep space 
increases (see Blum et al. 2015, for details). 

For example, the smallest aggregate size increases from 
2 mm to 6 mm for H 2 O ice activity at 1 AU if the thermophysi¬ 
cal model is used instead of the approximation. At 2 AU and for 
H 2 O ice sublimation, the smallest aggregate size changes from 
18 mm (approximation) to 174 mm (thermophysical model). 

The calculations shown in Fig. 1 were only performed for 
one aggregate size of 1 mm because the systematic simulations 
showed that starting from an aggregate radius of 15 pm, the vari¬ 
ation of the aggregate size does not considerably influence the 
resulting pressure at the ice-dust interface (see Fig. 2). 

With the knowledge of the gas pressure at the ice-dust inter¬ 
face, it is possible to derive the smallest dust-aggregate size that 
can be released from the surface by comparing the gas pressure 
with the tensile strength of the dust-aggregate layer. 

Strictly speaking, the assumption that all absorbed energy is 
used by the sublimation of the ice gives us the highest effective 
sublimation rate (i.e., highest mass loss per second) and not pres¬ 
sure. However, numerical tests agree with the values retrieved 
from Eq. 3 (see Fig. 1). Considering the thermal emission in 
the energy budget leads to a reduction of the effective loss mass. 
The role of this dissipation channel was discussed by Blum et al. 
(2015), where it was shown that this effect is weak close to the 
Sun and increases with heliocentric distance. Because the tensile 
strength of the covering aggregate layer is a decreasing function 
of aggregate size (see Eq. 1), the larger grains can be lifted eas¬ 
ily and their size is only restricted by the cometary gravity. In 
addition, the gas drag is stronger for larger aggregates because 
their cross section is larger. Thus, the application of Eq. 2 to¬ 
gether with Eq. 3 yields the smallest dust-aggregate size that can 
be released by the gas pressure at the ice-dust interface. For the 
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presented cases shown in Fig. 1, the thermophysical model pro¬ 
vides lower values for gas pressure and, hence, higher values for 
the smallest dust-aggregate size, because it takes into account 
other channels of the incoming energy dissipation. Nonetheless, 
for H 2 O ice, CO 2 ice, and CO ice, the approximation of the high¬ 
est gas pressure (Eq. 3) seems reasonable (see Fig. 1). 

The dust aggregates are lifted up by the outward-directed 
gas-friction force, which can be approximately described by 
Fgas = Jis^ piRjr)^ (Blum et al. 2015). Obviously, to re¬ 
lease aggregates from the cometary surface, the gas drag force 
has to exceed the inward-directed gravitational force, Fgrav - 
G (m M / R^) {Rjr)^, that is, Fgas ^ Fgrav Here, R is the radius of 
the comet nucleus (for simplicity, we assumed a spherical shape 
for the nucleus), r is the distance between the dust aggregate and 
the center of the nucleus, G is the gravitational constant, m is 
the mass of the dust aggregate, and M is the mass of the nucleus. 
Neglecting cohesion, we calculated the largest aggregate that can 
be ejected from the nucleus, 


tnux ^ ^ ip * \ / 

i O TT Cj Pcomet Paggregate ^ 

Here, pcomet and paggregate are the mass densities of the comet and 
the dust aggregates. We used an aggregate density of 875 kg m“^ 
(volume-filling factor of the aggregates multiplied by the bulk 
density). Equation 4 is only valid under the crude assumption 
that the nucleus and the released dust aggregate are spherical in 
shape and that the gas release is constant over the surface of the 
comet nucleus. Dust aggregates larger than the estimated size in 
Eq. 4 cannot be ejected into space. 


3. Application of the model to comet 
67P/Churyumov-Gerasimenko 

In this section, we apply the model described in Sect. 2 to 
comet 67P/Churyumov-Gerasimenko to derive the size range of 
the dust aggregates able to escape from the nucleus by the out- 
gassing of icy materials. The lower limit of the dust-aggregate 


size is determined by the tensile strength of the material that has 
to be exceeded by the gas pressure at the ice-dust interface (Eqs. 
2 and 3). The upper limit is given by the fact that the released 
aggregates have to escape the gravity field of the nucleus by the 
gas drag (outward-directed gas-friction force; see Eq. 4). 

The calculations were performed for H 2 O ice, CO 2 ice and 
CO ice. Each icy material was treated by a one-species model, 
which means that all available energy was consumed by the sub¬ 
limation of this particular icy species. To derive the largest dust- 
aggregate size, a radius of 2 km (Lowry et al. 2012) and a density 
of 470kgm“^ (Sierks et al. 2015) were used for the nucleus of 
comet 67P/Churyumov-Gerasimenko. Furthermore, we used a 
Bond albedo of 0.04, which is a typical value for comets (the ge¬ 
ometric albedo of comet 67P/Churyumov-Gerasimenko is 0.05; 
Sierks et al. 2015). 

Figure 3 shows the derived dust-aggregate sizes that can be 
released from the cometary surface and ejected into space by the 
sublimation of the volatile components as a function of helio¬ 
centric distance. The blue colors visualize at which heliocentric 
distances the dust activity can be expected with the outgassing- 
activity of the different ices. For comparison, the aphelion and 
perihelion of comet 67P/Chutyumov-Gerasimenko are denoted 
by the vertical dashed lines. 

As can be seen in Fig. 3 and as was reported in Blum et al. 
(2015), outgassing of H 2 O ice can only explain the dust activity 
within 2.5 AU. At larger heliocentric distances, the outgassing 
rate of H 2 O ice is not sufficient to release dust aggregates off the 
surface, hence, other more volatile materials must drive the dust 
activity at heliocentric distances larger than 2.5 AU. 

CO 2 ice and CO ice have a lower latent heat than H 2 O ice, 
therefore the saturation gas pressure of more volatile ices is 
higher than the corresponding H 2 O ice pressure if the energy 
available for the sublimation is the same. Thus, at a fixed he¬ 
liocentric distance, the activity of CO 2 ice or CO ice can release 
smaller particles than the outgassing of H 2 O ice (see the differ¬ 
ent blue areas in Fig. 3 for details). Additionally, the higher 
gas pressure caused by the sublimation of CO 2 ice and CO ice 
allows larger dust aggregates to escape from the gravity field of 
the nucleus. Thus, the size range of released particles expands 
in both directions. 

The outgassing of CO 2 ice and CO ice can explain the dust 
activity at a heliocentric distance of up to ~ 4 AU and ~ 5.5 AU, 
respectively. Thus, the gas-driven dust activity of comets is in 
general possible at large heliocentric distances. However, the 
size range of the released dust aggregates becomes smaller with 
increasing distance to the Sun and ends with a dust-aggregate 
size of ~ 1 dm close to the heliocentric distance at which the 
gas-driven dust activity vanishes. 

To explain a sustained activity at larger heliocentric distances 
(> 5.5 AU), other more volatile materials such as N 2 ice (latent 
heat of sublimation': A = 2.2 x lO^Jkg ') or O 2 ice (latent 
heat of sublimation: A = 2.3 x 10^ Jkg ') must be present close 
to the surface of the nucleus (for comparison, the latent heat of 
sublimation of CO ice is A = 3.0 x 10^ Jkg *). 

Figure 3 also shows that the size range of the dust aggregates 
that can be released off the nucleus widens when the comet ap¬ 
proaches the Sun, because the pressure at the ice-dust interface 
increases and thereby smaller dust aggregates can be released 
from the surface and larger dust aggregates are capable to leave 
the gravitational field of the nucleus. 

Recent observations performed by the Rosetta spacecraft 
(between 3.4 AU and 3.7 AU) in the inner coma of comet 

' http://encyclopedia.airliquide.com/Encyclopedia.asp?GasID=5 
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Figure 3. Radii of the dust aggregates whose adhesion can be overcome by the gas pressure (see Eqs. 2 and 3) and which are able to escape from 
the gravity field of the nucleus (see Eq. 4) of comet 67P/Churyumov-Gerasimenko as a function of heliocentric distance. The different blue colors 
denote the size range of the dust aggregates that are able to escape from the comet in case of outgassing of H 2 O ice, CO 2 ice, and CO ice. The 
calculations are performed for a radius of the nucleus of 2km (Lowry et al. 2012), an albedo of 0.04, and a density of 470kgm“^ (Sierks et al. 
2015). For comparison, the three particle populations (indicated by the numbers 1 to 3; see main text for details) measured by the Grain Impact 
Analyzer and Dust Accumulator (GIADA) and the OSIRIS Narrow Angle Camera (NAC) are shown by the dashed-dotted box and the dotted 
boxes (Rotundi et al. 2015). The smallest particle size of population 3 is 28 pm and is not shown in this figure). To derive the smallest aggregate 
size, we used a constant ice temperature that was determined by the exact thermophysical model at 1 AU for H 2 O ice (208.4 K; see caption of Fig. 
1), at 2 AU for CO 2 ice (117.2 K), and at 3 AU for CO ice (33.9 K; these positions are indicated by the stars). 


67P/Churyutnov-Gerasitnenko revealed different populations of 
particles close to the nucleus (Rotundi et al. 2015). Relatively 
large grains (between 2 cm and 1 m in radius; see the dashed- 
dotted box indicated with number 1 in Fig. 3) in bound or¬ 
bits have been detected by the OSIRIS Narrow Angle Camera 
(NAC). Figure 3 shows that the outgassing of CO ice can explain 
the presence of dust aggregates with radii in excess of 2 cm. 
However, larger aggregates (> 3 dm) cannot escape from the 
gravity field of the comet at 3.4 AU or 3.7 AU. The larger ag¬ 
gregates can have been released closer to the Sun (during former 
perihelion passages). For example, the ejection of meter-sized 
dust aggregates is only possible close to perihelion. At peri¬ 
helion of comet 67P/Churyumov-Gerasimenko, the derived size 
range fits the size range observed between 3.4 AU and 3.7 AU 
relatively well. 

The second population of particles detected by the OSIRIS 
NAC possesses radii between 1.3 mm and 1.3 cm (see the dotted 


box indicated with number 2 in Fig. 3). Additionally, the third 
population measured by the Grain Impact Analyzer and Dust Ac¬ 
cumulator (GIADA) consist of smaller particles with radii rang¬ 
ing from 28/tm to 262pm (see the dotted box indicated with 
number 3 in Fig. 3; the lower end of this particle population 
is not shown in the plot). The gap between these two popula¬ 
tions (2 and 3) is most probably caused by the detection limits of 
the two different instruments used to measure the particle size. 
This would imply that only parts of the particle size distribution 
in the inner coma of comet 67P/Churyumov-Gerasimenko have 
been investigated with these two instruments. 

The presented model of the gas-driven dust activity cannot 
explain the release of the particles found in populations 2 and 
3. The release of small dust aggregates (below the blue areas 
shown in Fig. 3) or of micrometer-sized monomer grains is not 
possible because of the relatively high tensile strength required 
to detach these particles from the surface (~ 1 kPa; see Blum 
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et al. 2014). Thus, the particles found in populations 2 and 3 
must be produced by disintegration of previously ejected larger 
dust aggregates. We recall that the calculations shown in Fig. 
3 are only valid for dust aggregates. For H 2 O ice aggregates or 
mixtures of dust and ice, the density of the aggregates decreases 
and the tensile strength between the aggregates most proabably 
increases, which slightly changes the result. However, if the 
Rosetta spacecraft were to detect icy particles, this would indi¬ 
cate that most probably CO ice or CO 2 ice drives the dust activity 
(the term dust includes H 2 O ice at this heliocentric distance) of 
comet 67P/Churyumov-Gerasimenko. 

The question remains which physical process is responsible 
for the production of the smaller dust-size fractions observed in 
comet 67P/Churyumov-Gerasimenko and in the comae and tails 
of other comets. Because of the estimated high internal tensile 
strengths of the released dust aggregates of > 100 Pa (Blum et al. 
2006), the stress required for the disintegration or fragmentation 
must exceed the tensile strength. Of the three obvious processes, 
(1) collisions among the dust aggregates, (2) electric charging of 
the dust aggregates, and (3) spin-up and centrifugal mass loss of 
the dust aggregates, respectively, collisions seem to be too infre¬ 
quent, and a typical electrostatic potential in the interplanetary 
space of 5 V is far too low to reach the required stresses. Thus, 
only rotational spin-up can be considered as a candidate for mass 
loss of the dust aggregates on their departure from the nucleus. 
Assuming spherical homogeneous dust aggregates of radius s 
and mass density Paggregate, spinning at a surface velocity v, the 
centrifugal stress along the equatorial region of their surface is 


with m and A are the mass and cross section of the part of the 
aggregate considered to be fragmented off. Assuming the latter 
to be spherical with radius io and a mass density po, we derive 



with Po = Paggregate ~ 1000 kg m For a spinning velocity of 

V ss 1 ms“' and sq ~ s, this stress approximately equals the typ¬ 
ical internal strength of the aggregates. For much higher spin¬ 
ning velocities of v ss 100ms“\ micrometer-sized monomer 
grains with tensile strengths of 10^ Pa (Heim et al. 1999) can 
be released if so/s > 10“^, that is, if s > 1mm. Future stud¬ 
ies should therefore concentrate on the gas-friction dynamics of 
the released dust aggregates with respect to the highest spin- 
up velocities. First empirical data on the rotation rate of emit¬ 
ted dust particles for comet 67P/Churyumov-Gerasimenko have 
been presented by Fulle, M., et al. (2015), who analyzed the 
light curves of more than 1,000 dust grains. They found that 
the most likely spin frequency of the particles is below 0.15 
Hz, with only a few aggregates rotating at a rate higher than 
1 Hz. They concluded that their observed particles had sizes 
in the millimeter range so that the highest rotation velocity was 

V ss 0.01 ms '. Thus, following Eq. 6, the highest centrifugal 
stress is (Tmax ~ 0.1 Pa, much too low to explain the formation 
of small dust aggregates by centrifugal splitting. 

4. Summary and discussion 

The main findings of this work can be summarized as follows 

- We developed a simplified analytical model to derive the size 
of the dust aggregates that can be released from the cometary 


surface into space by the sublimation of volatile materials 
covered by a nonvolatile surface layer. This model takes the 
cohesion of the material into account, which yields, together 
with the gas pressure at the ice-dust interface, the smallest 
size of the dust aggregates. Additionally, the gravitational 
force in combination with the gas friction force determines 
the largest size of dust aggregates that are able to escape from 
the gravity field of the comet. Our model is a one-species 
model, which means that all incoming energy is consumed 
by the sublimation of only one volatile component (in reality, 
comets are multicomponent mixtures with energy partition¬ 
ing among the different volatile species). Furthermore, this 
model is based on the assumption that comets have formed 
by gravitational collapse of dust agglomerates in the early 
solar system (Johansen et al. 2007; Skorov & Blum 2012; 
Blum et al. 2014), which implies that the surface of the nu¬ 
cleus consists of dust aggregates (the aggregates themselves 
are composed of micrometer-sized grains). We showed that 
the simplified analytical model agrees well with numerical 
calculations using a sophisticated thermophysical model (see 
Fig. 1 in Sect. 2). 

- The simplified model was then applied to comet 
67P/Churyumov-Gerasimenko, taking into account three 
different volatile materials: H 2 O ice, CO 2 ice, and CO ice 
(see Fig. 3 in Sect. 3). H 2 O ice can explain the gas-driven 
dust activity of comet 67P/Churyumov-Gerasimenko within 
a heliocentric distance of 2.5 AU. For CO 2 ice and CO ice, 
the gas-driven dust activity is possible up to ~ 4AU and 
~ 5.5 AU, respectively. Thus, outgassing of H 2 O ice cannot 
drive the observed dust activity between 3.4 AU and 3.7 AU 
(position of comet 67P/Churyumov-Gerasimenko during 
the first dust size measurements in the inner coma; Rotundi 
et al. 2015). 

- The size range of the dust aggregates able to escape from the 
nucleus into space widens when the comet approaches the 
Sun and narrows with increasing heliocentric distance (see 
Fig. 3 in Sect. 3). Independent of the volatile material, 
the size range always ends at larger heliocentric distances 
with an aggregate radius of ~ 1 dm, which is also the ob¬ 
served size of cometary trails (see, e.g., Reach et al. 2000; 
Harmon et al. 2004; Nolan et al. 2006; Ishiguro 2008; Trigo- 
Rodriguez et al. 2010). The size range of the dust aggre¬ 
gates widens when the comet approaches to the Sun because 
the tensile strength of the dust aggregates decreases with in¬ 
creasing aggregate size. Thus, with higher pressure avail¬ 
able at the ice-dust interface closer to the Sun, smaller dust 
aggregates can be released from the surface. Furthermore, 
the pressure increase enables larger dust aggregates to es¬ 
cape against the gravity of the nucleus into space. This can 
explain the presence of the aggregates larger than ~ 3 dm 
if these aggregates have been released during former peri¬ 
helion passages. Recently, Soja, R. H., et al. (2015) found 
that the trail of comet 67P/Churyumov-Gerasimenko is dom¬ 
inated by millimeter-sized particles. Our model can explain 
the release of these particles, but only with the activity of CO 
ice close to perihelion. 

- The model can explain the presence of bound grains (with 
radii between 2 cm and 1 m ; Rotundi et al. 2015) close to the 
nucleus of comet 67P/Churyumov-Gerasimenko (dashed- 
dotted box in Fig. 3). However, the smaller particles detected 
by GIADA and OSIRIS NAC onboard the Rosetta spacecraft 
at heliocentric distances between 3.4 AU and 3.7 AU (dotted 
boxes in Fig. 3; Rotundi et al. 2015) and between 2.4 AU and 
3.4 AU (Della Corte, V., et al. 2015) cannot be released by 
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the gas pressure at these heliocentric distances. These par¬ 
ticles must be the product of a disintegration of larger dust 
aggregates in the coma. Furthermore, our model roughly 
agrees with the findings by Fulle, M., et al. (2015), show¬ 
ing that particles with radii of a few millimeters are a good 
model for the OSIRIS observations of free-floating particles 
in the coma. 

- The release of micrometer-sized monomer grains is not pos¬ 
sible because of the relatively high tensile strength required 
to detach the grains from the surface (~ 1 kPa; see Blum 
et al. 2014). This high gas pressure cannot be reached at the 
ice-dust interface under normal conditions. 

The results we discussed can also be used to model the gas- 
driven dust activity of other comets. The application of our 
simplified model to another comet would shift the largest dust- 
aggregate size (upper limitation of the blue areas in Fig. 3) and 
thereby the heliocentric distance at which the dust activity starts. 
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